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ABSTRACT 

The  work  done  in  this  investigation  is  part  of  a  continuing  program 
designed  to  develop  probabilistic  strength  and  life  models  for  filament 
composite  materials.   In  this  particular  experiment,  groups  of  single  fibers 
from  two  chemically  identical  different  production  spools  of  graphite  were 
loaded  to  identical  sustained  tensile  loads  in  order  to  determine  their  life 
statistics  while  under  stress.   While  previous  work  has  produced  models 
of  composite  reliability  as  a  function  of  strength,  this  work  develops  a 
model  which  will  be  necessary  to  predict  the  reliability  of  composite  fibers 
as  a  function  of  time.    Parametric  and  non-parametric  methods  were  used 
to  quantify  the  relationship  between  fiber  reliability  as  a  function  of  stress 
and  as  a  function  of  time. 
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I.    INTRODUCTION 

While  still  in  its  infant  stages,  the  use  of  composite  materials  in 
manufacturing  is  steadily  increasing,  particularly  in  high  reliability,  low 
accessibility  environments  such  as  aircraft  fuselages  and  satellite  fuel  tanks. 
High  strength  to  weight  ratios,  combined  with  minimal  maintenance 
requirements  and  extreme  resistance  to  corrosion  make  composites  ideally 
suited  for  structures  previously  built  exclusively  of  steel  or  aluminum. 
However,  because  composite  usage  has  such  a  short  history,  there  is  no 
significant  data  base  which  would  allow  accurate  reliability  predictions. 

Through  an  ongoing  program  of  research  and  experimentation,  The 
Advanced  Composites  Laboratory  at  the  Naval  Postgraduate  School  at 
Monterey,  CA,  is  defining  a  probabilistic  model  which  is  based  on  the 
failure  mechanism  of  composites  in  tension.  Small  scale  designs  which 
exploit  fiber  strength  under  tension,  have  accurately  verified  the  Local 
Load  Sharing  model.   More  recent  work  is  concentrating  on  the  application 
of  this  model  to  large  scale  structures.   With  this  in  mind,  it  is  important  to 
understand  that  a  structure  must  not  only  be  predictable  (reliable)  in 
strength,  but  also  in  life.   Because  of  their  non-homogeneity,  the  study  of 
composite  materials  must  be  further  broken  down  to  the  study  of  not  only 
individual  constituents,  but  also  of  their  inter-relationships  with  each  other. 
For  our  purposes  here,  all  work  is  done  exclusively  with  graphite/  epoxy 
composite.;  therefore,  a  proper  understanding  of  the  fiber,  the  matrix,  and 
the  interaction  between  them,  is  essential  to  formulating  an  accurate  model 
of  composite  reliability. 


The  primary  objective  of  this  project  was  to  examine  a  large  sample 
space  of  graphite  fibers  under  different  load  conditions,  and  formulate  a 
model  of  their  reliability  in  life.   The  fiber  life  results  would  then  be 
applied  to  composite  life  experiments  now  on-going  designed  to  predict  the 
life  of  composite  structures.   With  data  from  these  two  projects  in  hand, 
one  could  then  ostensibly  isolate  the  solution  to  the  fiber/  matrix  interface 
problem  and  produce  a  realistic  life  prediction  model. 


II.  BACKGROUND 

A.     STRENGTH  DETERMINATION  OF  FIBER  REINFORCED 
COMPOSITES 

Fiber  Reinforced  Composite  materials  generally  consist  of  high 

strength  fibers  bonded  together  in  a  ductile  matrix.    In  these  materials  the 

fibers  carry  essentially  all  of  the  stresses  imposed  on  the  structure  while 

the  matrix  serves  as  the  primary  vehicle  for  transferring  load  between 

adjacent  fibers,  known  as  the  Local  Load  Sharing  model.   In  this  model, 

the  local  load  carried  by  a  bundle  of  fibers  is  carried  equally  until  the 

weak  fibers  begin  to  rupture.   When  a  fiber  breaks,  its  load  is  transferred 

to  the  adjacent  fibers  through  the  matrix.   Rosen  has  modeled  the  load 

transfer  by  the  existence  of  an  ineffective  length  at  the  end  of  each 

(failed)  fiber  [Ref.   1].  In  this  region,  the  longitudinal  stress  (a)  diminishes 

to  zero  while  shear  stresses  (x)  rise  to  a  maximum.   As  these  ineffective 

lengths  accumulate,  loads  continuously  transfer  from  the  weak  (failed) 

fibers  to  the  strong  ones.   If  enough  fiber  failure  sites  are  in  close 

proximity  to  each  other,  the  matrix  begins  to  crack  or  de-bond  from  the 

fibers.     When  this  occurs,  the  matrix  can  no  longer  transfer  the  load  and 

the  sample  fails  catastrophically.   Thus  composite  strength  can  be 

modeled  as  a  function  of  fiber  strength,  the  ability  of  the  matrix  to  transfer 

load,  and  the  interface  between  the  fiber  and  matrix.   The  mathematical 

model  which  predicts  the  probability  of  failure  of  the  composite  (in 

tension)  from  the  statistics  of  fiber  strength  and  the  effective  length  has 

been  established  by  Harlow  and  Phoenix  [Ref  2  and  3]. 


These  three  quanta  are  recognized  as  building  blocks  in  composite 
strength  determination.  Matrix  optimization  is  guided  by  this  model.  From 
a  given  assortment  of  matrix  materials,  one  may  determine  which  has  the 
most  suitable  characteristics  for  stiffness  and  transferring  shear  stresses. 
Fiber  strength  statistics  are  determined  through  loading  to  destruction 
measurement.   However,  it  is  not  feasible  to  directly  measure  and  quantify 
the  molecular  level  interface  region.   This  can  be  accomplished  by 
measuring  the  strength  of  a  composite  of  known  fiber  and  matrix,   the 
fiber  strength  statistics  and  composite  strength  statistics  can  then  be  used 
simultaneously  to  solve  for  the  ineffective  length  parameter. 

B.     LIFE  DETERMINATION  OF  FIBER  REINFORCED 
COMPOSITES 

Decades  of  experience  with  homogeneous  materials  support  the 

correlation  between  strength  and  life:  that  is,  the  larger  the  factor  of 

safety  for  a  given  structure,  the  longer  it  can  be  expected  to  last. 

However,  Coleman's  theory  of  time  dependence  in  the  mechanical 

breakdown  of  fibers  [Ref.  4]   implies  that  the  strength/life  relation  in  load 

bearing  fibers  may  not  be  so  straightforward.   Using  the  "weakest  link 

hypothesis"  first  developed  by  Pierce  in  1926,  Coleman  concluded  that 

fiber  life  was  instead  a  function  of  the  statistical  distribution  of  these 

"weak  links"  in  time  dependent  strength.   The  time  dependent  strength  of 

the  fiber,  termed  breaking  kinetics,  is  the  extreme  value  of  fiber  life  under 

a  given  load  history.  This  is  fundamental  to  the  modeling  of  composite  life. 


Only  when  fiber  life  statistics  are  available  can  the  model  be  applied  to 
composite  structures.   Hence,  the  purpose  for  this  work. 

The  theory  proposed  is,  the  strength  of  a  specific  fiber  can  be  related 
to  the  life  of  that  same  fiber.  The  dichotomy  of  course  is,  once  a  fiber  has 
been  broken  in  the  strength  test,  another  molecularly  identical  sample  is 
not  available  for  the  life  test.  The  alternative  is  to  obtain  a  set  of 
statistically  identical    fibers  and  test  a  portion  of  them  in  strength  and 
another  portion  in  life. 

This  is  being  accomplished  by  using  two  different  spools  of  AS-4 
graphite  fiber,  AS4-008  and  AS4-019.   A  statistically  meaningful  number 
of  samples  have  been  taken  from  these  two  spools  to  estimate  the 
respective  parameters  of  the  Weibull  Distribution.   The  estimated 
parameters  of  these  two  different  spools  are  different,  suggesting  the 
intrinsic  strength  of  these  two  spools  are  from  two  different  populations. 


III.     FILAMENT  LIFE  TESTING  METHOD 

A.  TEST  SPECIMENS 

Test  samples  were  single  strands  of  Hercules  Magnamite  AS4  cut  to 
two  inch  gauge  lengths.   Each  sample  consisted  of  one  fiber  filament 
mounted  on  a  paper  cardboard  with  a  cutout  for  handling.   The  cardboard 
mount  was  severed  with  a  hot  wire  prior  to  load  application.   Specific 
handling  and  mounting  procedures  are  described  in  Appendix  A  and  in 
Reference  5.    Once  mounted,  the  fiber  diameters  were  measured  using  a 
laser  diffraction  device  also  following  the  procedures  described  in 
Reference  5. 

B.  TEST  EQUIPMENT  DESCRIPTION 

The  single  filament  samples  were  arranged  in  groups  of  64  at  four 
different  sustained  load  levels  by  dead  weights.   The  elapsed  time 
between  the  startup  of  load  and  ending  of  stress-rupture  at  each  station 
was  triggered  by  an  infra-red  light  switch.   Test  equipment  is  described  in 
Appendix  B.    A  polymethyl  methacrylate  (PMMA)  cover  was  used  to 
protect  the  samples  from  ambient  air  disturbances  and  ultraviolet  (UV) 
radiation.   It  is  known  that  UV  light  degrades  the  strength  of  some  fibers; 
therefore,  PMMA  was  chosen  as  a  shield  and  is  analyzed  in  Appendix  C. 
Once  loaded,  the  fibers  were  continuously  monitored  by  a  dedicated  PC 
(personal)  class  computer.    Software  programs  are  reproduced  in 
Appendix  D. 


C.     TESTING  PROCEDURE 

1.  Strength  Test 

Fiber  strength  data  were  acquired  previously  and  are  covered  in 
References  5  through  7. 

2.  Life  Test 

During  the  diameter  measurement  process,  every  fifth  fiber  was 
tested  to  failure  in  strength  using  a  calibrated  load  cell.  This  provided  the 
link  between  the  two  groups  needed  to  verify  that  the  samples  tested  in 
life  were  statistically  identical  to  those  tested  previously  in  strength.   The 
actual  loading  of  the  samples  on  the  board  was  carefully  controlled  and, 
as  much  as  possible,  were  all  done  by  the  same  operator.   The  loading 
procedure  is  also  described  in  Appendix  A. 

There  were  four  levels  of  load  set  up  on  the  board,  14.17  g. 
(Level  1),  13.15  g.  (Level  2),  and  12.20  g.  (Level  3)  and  10.96  g.  (Level 
4).   Each  load  level  had  124  stations,  64  for  each  spool.  This  allowed  a 
total  of  512  life  samples  to  be  measured  at  sustained  loads.   Level  4  was 
planned  but  postponed.    Instead  those  stations  were  also  loaded  to  14.17 
g.  to  gain  additional  measurements  at  that  level.   The  sample  board 
alternated  rows  with  008  and  019  fiber  to  give  4  rows  of  each.   After 
each  fiber  was  loaded,  the  date  and  time  of  the  load  were  recorded  in  the 
lab  notebook  and  also  in  a  software  spreadsheet.   If  a  fiber  failed  on 
loading,  its  failure  stress  was  recorded.    After  loading,  each  station  was 
actively  monitored  by  an  AT  computer  which  recorded  stress  failures  as 
they  occurred.   Exact  failure  information  was  periodically  retrieved  and 
recorded  in  the  notebook  and  spreadsheet.    The  spreadsheet 


automatically  computed  realized  data  for  stress  ruptures,  early  censor 
(Type  Iv)  stress  ruptures,  fail  on  load  stresses,  and  existing  lifetimes  for 
each  fiber  type  at  each  load  level.   It  would  then  arrange  the  desired  data 
in  a  form  easily  transported  to  other  programs  for  analysis  or  presentation. 


IV.     RESULTS 

A.     TEST  DATA 
1.     Overview 

The  data  presented  herein  will  be  used  to  identify  a  suitable 
strength-life  model  for  fiber  filaments.  A  probabilistic  formulation  of  the 
strength-life  model  is  a  joint  distribution  of  the  random  variables,  stress 
(applied),  a,and  elapsed  time,  t,  denoted  f(o,t).  The  parameters  of  this 
joint  distribution  are  determined  by  three  classes  of  realized  data.   The 
first  is  data  realized  in  strength  over  negligible  time.   The  second  is  data 
realized  in  life  for  a  given  stress  history,  in  this  experiment,  constant, 
sustained  stress.   The  third  is  data  realized  in  strength  given  a  survival  in 
time  history.   The  first  two  cases,  that  is,  strength  and  life  data,  will  be 
presented  in  this  investigation.  To  examine  variability  in  production,  two 
spools  of  AS-4  graphite  fibers  were  treated  as  two  distinct  populations, 
designated  as  spool  008  and  019.  The  division  of  fibers  from  these  two 
spools  is  traced  out  in  Figures  1  and  2.  The  entire  population  of  either 
spool  was  denoted  as  {D}.   This  population  is  divided  into  strength  test 
samples  yielding  data,  {Ds },  and  life  test  samples  yielding  data,  {D^  }.   The 

respective  sub-division  of  these  two  data  sets  is  described  in  the  following 
two  sections. 


2.     Strength  Data 

The  strength  data  set  comes  from  work  done  previously  using  the 
same  008  and  019  fiber  populations  [Refs.   5  through  7].   The  sources  of 
this  data  are  summarized  in  Table  1. 

TABLE  1.    SOURCES  OF  FIBER  STRENGTH  DATA 
DATA  NAME  REF.  NO. 


SET 

{Ds1} 

Dr.  Edward  M.  Wu 

6 

1      S  1  ' 

Mr.  Glen  Nypiak 

(Ds2} 

Lt.  David  Bell,  USN 

5 

{Ds3} 

LCDR  Carl  Englebert, 

USN 

7 

{Ds4  }  Mr.  Jim  Nageotte 


3.     Life  Data 

Within  the  life  data  set,  {D-^  },  not  all  of  the  samples  were  realized 

in  life.    Some  of  the  samples  were  tested  to  assure  the  life  samples  were 
statistically  the  same  as  the  strength  samples  and  are  denoted  as  {D-^s} 

(Figures  1  and  2).  Some  samples  failed  before  they  reached  the  target 
sustained  load  level.  Thus,  this  was  part  of  the  life  test  subset  but  was 
realized  in  strength  and  is  denoted  as  {D^s},  {D-^sh  and  {D-^}.  Based 

on  kinetic  theory  [Ref.  2],  it  is  generally  believed  that  the  load  supported 
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by  a  fiber  increases  as  the  loading  rate  increases;  therefore,  failures  in 
strength  would  be  hypothetically  eliminated  if  the  fibers  could  be  loaded 
instantaneously.  While  this  not  physically  possible,  the  data  from  samples 
which  failed  on  loading  can  be  accounted  for  by  treating  the  life  samples 
conditionally  to  strength  greater  than  the  sustained  load.    Data  obtained 
from  this  experiment  is  summarized  in  Tables  2  through  4.   Ruptures 
occurring  on  loading,  prior  to  the  desired  load  are  identified  as  "Fail  on 
Load".    Stress  ruptures  occurring  in  time  but  caused  by  other  than 
intrinsic  stimulus  internal  to  the  sample  are  "Type  Iv"  censored.    This  data 
is  statistically  treated  as  greater  than  or  equal  to  the  realized  life.   The 
censorings  are  chance  events  as  denoted  by  the  variable  time  at  censor. 
Stress  ruptures  in  time  are  "Exact"  and  fibers  still  under  load  at  the  end  of 
the  experiment  are  "Type  I"  censored.    This  data  is  statistically  treated  as 
greater  than  or  equal  to  the  time  the  data  was  interpreted. 

TABLE  2.  RESULTS  OF  LOAD  LEVEL  1 


Type  of  Failure 

AS4-008 

AS4-019 

Fail  on  Load  { D^  j  s } 

53 

57 

Type  I v  Censor  { D-j-  1v } 

9 

7 

Exact  {Df-u  } 

15 

20 

Type  I  Censor  {D^c} 

116 

110 

11 


TABLE  3.  RESULTS  OF  LOAD  LEVEL  2 


Type  of  Failure 

AS4-008 

AS4-019 

Fail  on  Load  {D-^s } 

26 

13 

Type       Iv       Censor 

4 

1 

(Dt2v) 

Exact  {D^2t } 

9 

6 

Type  I  Censor  {D^c  } 

62 

63 

TABLE  4.  RESULTS  OF  LOAD  LEVEL  3 


Type  of  Failure 

AS4-008 

AS4-019 

Fail  on  Load  {D^  } 

10 

11 

Type  Iv  Censor 

1 

3 

{Dt3v} 

Exact  {Dt3t} 

3 

1 

Type  I  Censor  {D^^  } 

62 

63 

B.   DATA  ANALYSIS 

1.     Analysis  as  a  Function  of  Stress 
a.      Overview 

Before  proceeding  with  an  analysis  of  life  based  on  strength, 
it  must  first  be  confirmed  that  the  two  sets  of  data,  {Ds }  and  {D-^  },  are 

indeed  statistically  from  the  same  population.   As  can  be  seen  in  Figures  1 
and  2,  each  spool  has  yielded  two  sets  of  samples.   The  first  set,  {Ds }, 
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was  tested  entirely  in  strength.   Within  the  second  set,  {D-j.  },  every  fifth 
fiber  was  tested  in  strength  in  order  to  link  the  two  sample  sets  {Ds }  and 

{Dtl- 

b.      N on- Parametric   Interpretation 

Non-parametric  data  analysis  was  conducted  by  plotting  the 
cumulative  failures  in  a  weakest  link  (Weibull  Probability)  coordinate. 
This  required  the  data  to  be  ordered  and  a  rank  assigned.   For  this 
experiment,  Expected  Rank  was  used  and  is  defined:  the  ith  realization  of 
the  ordered  Xj  is 

— —     for  i  =  1.2,3.. .N 
N  +  l 

and  Xj  is  the  ordered  data  (from  weakest  to  strongest)  in  a  set  of  N 
samples.   An  ambiguity  arises  when  one  considers  the  inclusion  of 
censored  data.    Because  the  point  is  censored,  its  exact  strength,  and 
therefore  its  exact  rank,  cannot  be  precisely  determined.    Furthermore, 
because  the  data  is  ranked,  this  uncertainty  perturbs  the  rank  of  each 
subsequent  data  point  greater  than  the  censored  data.    For  this  reason, 
censored  data  was  omitted  from  the  graphical  representation.    However, 
censored  data  did  not  affect  the  Likelihood  calculations  and  was  therefore 
included  in  the  parametric  interpretation. 

Once  the  data  has  been  ranked,  it  is  linearized  by  using  the 
weakest  link  transformation.   This  method  linearizes  the  ranked  data  and 
allows  a  visual  observation  of  the  Tower  tail  characteristics.   The  weakest 
link  transformation  is  defined  as: 

F*(o)  =  ln(-ln(l-F(o))) 
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where  F(o)  is  the  empirical  rank  of  the  data  in  question.   Recalling  that 
each  fiber  is  modeled  as  a  chain  of  independent  links,  the  strength  of  a 
given  length  of  fiber  must  be  the  same  as  the  strength  of  the  weakest  link 
in  that  chain,  hence  F*(o)  can  be  further  defined  as  the  probability  of 
failure  of  the  weakest  link  of  a  given  gauge  length  of  fiber. 

Strength  data,  {Ds },  from  these  two  spools  (008  and  019)  is 

plotted  in  Figures  3  and  4.  This  set  of  fibers  had  previously  been 
confirmed  to  have  come  from  the  same  set  of  like  materials  and  gauge 
lengths  (Ref.   9).   Figures  5  and  6  depict  the  F*  plot  of  the  set  D^s  overlaid 
with  the  set  Ds .   The  close  proximity  of  the  two  sets  of  data  indicates  they 
too  are  of  the  same  fiber  and  gauge  length.   A  cursory  run-test  on  the 
ranked  set  of  merged  data  shown  in  Figures  7  and  8  shows  acceptable 
levels  of  statistical  clumping  and  completes  the  non-parametric 
identicalness  test  of  the  two  sets. 

c.      Parametric   Interpretation 

A  parametric  analysis  of  the  same  data  was  also  conducted 
by  using  the  Weibull  distribution  for  both  strength  and  life: 


F(x)  =  1-exp 


(  r' 

x 


IP 

where  a  and  (3  are  the  shape  and  location  parameters  respectively. 
Given  a  data  set,  these  parameters  may  be  obtained  by  using  the 
Maximum  Likelihood  Estimator  (MLE)  derived  in  Appendices  E  and  F. 
The  attribute  of  this  method  is  there  is  no  requirement  to  order  the  data 
which  eliminates  the  process  of  ranking  which  allows  the  inclusion  of  any 
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censored  data  points.    Also  associated  with  the  MLE  is  the  concept  of 
confidence  interval.   The  confidence  region  of  the  parameters  a  and  p 
can  be  obtained  from  the  likelihood  contour  plots.  The  confidence  region 
is  the  confidence,  stated  as  a  percentage,  that  a  selected  a  and  p  are  the 
actual  a  and  p.     Parameters  chosen  at  a  very  low  confidence  interval 
have  a  low  probability  of  matching  the  actual  parameters,  but  when  they 
do,  they  are  highly  accurate.   A  very  high  confidence  interval,  while 
allowing  more  room  for  error,  is  not  very  precise  and  is  of  little  practical 
value. 

Figures  9  and  10  show  a  three  dimensional  representation  of 
the  Likelihood  for  the  Ds  data  set  along  with  various  confidence  region 

contours.   The  5  %  and  1  %  confidence  region  contours  are  shown  in 
Figures  11  and  12.   The  shape  parameter,  oc,  and  location  parameter,  p, 
for  this  data  set  can  be  determined  from  Figures  3  and  4  graphically  by 
the  slope  of  the  line  of  data  points  and  the  value  of  the  X  axis  at  F*  =  0 
respectively.    Alternately,  the  estimate  can  also  be  based  on  the  calculus 
MLE  derived  in  Appendix  E. 

In  order  to  show  that  {Ds }  is  identical  to  {D-^s},  the  likelihood 

ratio  can  be  used.  The  likelihood  ratio  is  graphically  represented  by  the 
volume  of  the  intersection  of  the  two  likelihoods  at  a  given  confidence 
interval,  divided  by  the  volume  of  the  same  confidence  interval  of  the 
accepted  baseline  data  set,  {Ds}  or,  symbolically 

V(D.)r|V(DJ 
V(D.) 
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:  V(DS)  and  V(D^S)  are  the  volumes  of  the  likelihood  contours 

within  the  desired  confidence  interval.   Figures  13  and  14  show  overlay 
plots  of  the  Likelihood  contours  of  Ns  and  N{s  at  5%  and  1  %  confidence 

intervals.    By  examining  these  figures,  it  can  be  determined  that  the 
likelihood  ratios  for  008  and  019  at  the  5  %  confidence  interval  are 
estimated  to  be  .75  and  .70  respectively,  which  is  deemed  to  be 
acceptable  for  this  confidence  interval. 

The  two  sets  of  data  can  now  be  merged  and  the  tested  in 
life  fibers  are  confirmed  to  be  of  the  same  population  as  the  tested  in 
strength  fibers.   As  a  last  point  of  confirmation,  Figures  15  and  16  overlay 
this  merged  data  set  {Ds+D^s}  with  the  original  set,{Ds}  showing  they 

are  essentially  identical,  as  do  the  contour  overlay  plots  shown  in  Figures 
17  and  18. 

A  distinction  can  now  be  made  concerning  those  fibers 
which  failed  on  load  ({D^1s},  {D^s  h  {D^  }).   Although  intended  to  be 

measured  in  time,  these  fibers  actually  failed  in  strength  prior  to  reaching 
the  target  loads,  c^ ,  G^  and  G3  respectively;  therefore,  this  set  of  data  can 

be  considered  as  realized  in  strength  to  be  included  in  the  merged  strength 
data  set  {Ds+D^s}.   The  Expected  Rank  of  the  exact  failed  in  strength 

data  was  presented  as  before  with  the  exception  that,  the  samples  which 
reached  the  target  load  in  life  test  were  considered  Type  I  strength 
censored  at  the  stress  rupture  load,   This  is  the  equivalent  to  stating  its 
strength  would  be  greater  than  the  stress  rupture  load.   As  with  all  non- 
parametric  analyses,  these  censored  points  were  not  presented  in  the  F* 
domain  because  of  the  ambiguity  in  rank.   By  treating  censored  data  in 
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this  manner,  the  graphical  representation  of  the  middle  and  upper  tails 
remains  unaffected,  while  the  resolution  of  the  lower  tail  is  improved.  Fail 
on  loading  data,  {D^},  {D^},  and  {D-^},  is  presented  in  Figures  19 
through  24  and  can  be  analyzed  in  the  same  manner  as  the  set  { D^  s }  with 
similar  conclusions  being  drawn.   Thus,  the  set  of  fibers  tested  in  time  can 
now  be  considered  statistically  identical  to  the  set  of  fibers  tested  in 
stress. 

The  finaJ  step  in  strength  analysis  was  to  analyze  the  entire 
population  of  strength  data  and  determine  a  shape  (a)  and  location  ((3) 
parameter  for  the  Weibull  Probability  of  Failure  Function  which  would  be 
necessary  to  determine  the  Likelihood  in  Life.   Results  of  non-parametric 
analysis  of  the  set  {D}  in  strength  are  shown  in  Figures  25  and  26. 
Because  of  the  data  censored  above  12.20,  13.15,  and  14.17  grams,  a 
value  for  a  and  (3  can  only  be  estimated  from  the  set  {Ds  +  D^s}; 

therefore,  a  parametric  analysis  in  the  form  of  the  Likelihood  contours  of 
the  set  {D}  is  shown  in  Figures  27  and  28. 
2.     Analysis  as  a  Function  of  Time 

Since  the  life  of  a  fiber  is  dependent  on  the  fiber's  surviving  the 
loading  process,  ruptures  during  loading  must  be  considered  in 
determining  the  Likelihood  of  the  fiber  in  life.  The  reliability  of  the  fiber  in 
strength  will  be  necessary  to  determine  its  reliability  in  life.   Therefore, 
equation  (1)  in  Appendix  E  must  be  modified  to  include  the  fibers 
reliability  in  strength: 
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R(o)  =  exp 
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and  the  Likelihood  in  the  time  domain  becomes: 
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where  oc  is  the  target  stress  and  k  is  the  number  of  left  censored  fibers. 

By  determining  estimates  of  (Xq  and  pa  from  Figure  27  for  Spool  008  and 
evaluating  the  Likelihood  contour  over  a  set  range  of  a^  and  p^ ,  the 

Likelihoods  for  each  of  the  three  load  levels  can  be  plotted  and  are  shown 
in  Figures  29  through  31  which  were  developed  using  an  estimate  of  5.0 
for  a<j  and  18.0  for  pa.  A  similar  procedure  was  followed  for  Spool  019. 

Although  not  investigated  during  this  analysis,  an  alternative 
method  of  calculating  the  Likelihood  would  be  to  evaluate  the  fail  on 
loading  data  as  censored  in  time  rather  than  stress.   In  order  to  do  this, 
one  must  first  assign  a  time,  tg,  at  which  the  sample  is  considered  to  have 

attained  the  target  stress.   All  samples  which  fail  on  loading  can  be 
considered  to  have  a  life  less  than  or  equal  to  tg,  or: 
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Under  this  interpretation,  the  Likelihood  can  be  calculated  using: 
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3.     Construction  of  the  Strength-Life  Model 

For  a  life  distribution  given  a  sustained  stress  level,  or  the  Weibull 

model  used  was: 


F(tlc  =  c.)  =  1  -  exp' 


<  t   ^ 


P 


VH,  ) 


Using  this  model,  the  Likelihood  contours  for  the  life  data  realized  up  to 
the  current  time  were  plotted  and  are  shown  in  Figures  29  through  31. The 
flatness  of  the  contours  is  the  result  of  an  extremely  broad  range  of  (3 
which  is  caused  by  the  elapsed  time  (at  the  time  of  data  interpretation) 
being  much  smaller  than  the  underlying  location  parameter  of  life,  (3. 
Therefore,  only  very  tentative  estimates  of  (3  are  possible  at  this  time. 
The  shape  and  location  parameters  are  visually  estimated  to  be: 
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Stress  Level,  Oj 

&, 

A 

14.17  g. 

.15 

1.0E+11 

13.15  g. 

.15 

1.0E+12 

12.20  g. 

.15 

5.0E+13 

The  estimated  life  at  different  levels  of  sustained  stress  is  shown 
in  Figure  32  evaluated  at  three  levels  of  probability  of  failure,  0.001,  0.01, 
and  0.1.   The  figure  can  be  considered  as  a  graphical  representation  of  the 
Strength-Life  Model  for  AS4  graphite  under  sustained  stress.   The 
expected  life  of  the  fiber,  evaluated  at  the  desired  probability  of  failure  at  a 
given  stress  level,  can  be  obtained  directly  from  the  model. 
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V.      CONCLUSIONS   AND   RECOMMENDATIONS 

The  Strength-Life  Model,  once  verified,  can  be  a  useful  and  important 
tool  for  predicting  fiber  reliability  using  strength  statistics.  In  design,  it 
allows  accurate  determination  of  maximum  design  stress  for  an  expected 
life.   In  maintenance  and  repair,  it  facilitates  the  calculation  of  the 
remaining  life  of  an  existing  structure.   It  also  predicts  stress  levels  at 
which  an  existing  structure  may  be  operated  that  can  increase  or 
decrease  its  remaining  life. 

The  following  items  are  recommended  for  follow-on  research:. 

1.  Application  and  verification  of  the  model  using  larger  scale 
composite  samples. 

2.  Mathematical  formulation  to  predict  life  based  solely  on  strength. 

3.  The  life  test  be  continued  to  improve  the  estimation  of  the 
parameters  of  the  model. 
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APPENDIX    A 
PROCEDURES. 

Fiber  Handling.   Fiber  handling  procedures  are  specific  and  are  in 
effect  continuously. 

1 .  Unless  required  for  the  immediate  task  at  hand,  all  fibers  and 
filaments  shall  be  stored  in  dark,  secure  containers. 

2.  Fibers  and  filaments  are  fragile.   Any  fiber  or  filament  that  is 
dropped  or  otherwise  subjected  to  extreme  stress,  shall  be  discarded. 

3.  Whenever  possible,  all  samples  should  be  handled  by  one  end 
only  so  as  not  to  exert  any  tension  on  the  sample. 

4.  Use  of  fluorescent  lights  shall  be  kept  to  a  minimum. 

Filament  Mounting.     Filaments  were  prepared  according  to  the 
following  guidelines: 

1.  Strip  off  a  random  length  of  fiber  from  the  spool,  sever,  and 
discard. 

2.  Cut  the  next  12-15  cm.  of  fiber  from  the  spool  and  tape  one  end 
to  the  glass  plate. 

3.  Pour  a  small  puddle  of  ethyl  alcohol  on  the  glass  plate  to  float  the 
free  end  of  the  bundle.  The  filaments  will  begin  to  float  apart.  Gently 
fan  out  the  filaments  floating  in  the  alcohol  and  let  stand  until  all  the 
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alcohol  has  evaporated. 

4.  Write  the  correct  serial  numbers  on  the  carriers  to  be  loaded. 

5.  Using  a  small  piece  of  transparent  tape  on  the  end  of  a  modelling 
knife,  lift  the  free  end  of  one  filament  and  slide  an  empty  carrier  under 
it 

6.  While  maintaining  filament  alignment  with  the  two  small  holes  in 
the  carrier,  lower  the  raised  filament  end  and  tape  it  to  the  glass 
plate. 

7.  Use  two  small  pieces  of  tape  to  secure  the  filament  to  the 
carrier. 

8.  After  the  filament  is  secured  between  the  two  pieces  of  tape, 
carefully  cut  the  filament  to  free  it  and  the  carrier  from  the  glass  plate 
and  place  the  filament/carrier  assembly  on  the  cold  curing  plate. 

9.  Repeat  steps  5-8  eight  more  times.    After  the  ninth  filament  has 
been  mounted,  discard  the  remains  of  the  bundle. 

NOTE 

In  order  to  maintain  standard  curing  times  and  temperatures,  do 
not  activate  the  heat  element  in  the  hot  plate  until  the  plate  is  full  and 
do  not  remove  the  samples  until  the  plate  has  been  turned  off  and  has 
cooled. 

10.  Activate  the  hot  plate  with  a  setting  of  200°  F. 

1 1 .  Thoroughly  mix  a  small  amount  of  epoxy  (Devcon  2  Ton). 

12.  Reduce  hot  plate  temperature  to  150°  F. 
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13.  Place  one  dot  of  epoxy  at  each  end  of  the  elongated  hole  in 
each  carrier  to  permanently  secure  both  ends  of  the  filaments  to  their 
respective  carriers.   Also,  place  a  dot  on  the  outboard  ends  of  the  two 
small  holes  to  re-enforce  them. 

14.  When  the  last  filament  has  been  epoxied,  turn  off  the  hot  plate 
and  let  cool. 

15.  Store  mounted  filaments  in  designated  containers. 

Filament   Loading. 

Loader  Start  up 

1 .  Turn  on  power  supplies  for  load  cell  and  elevator. 

2.  Turn  on  HP-85  and  attached  Data  Acquisition  Unit. 

3.  Insert  program  tape  into  HP-85  and  type  LOAD  "QUIK4". 

4.  Type  RUN. 

5.  Select  "k2"  (INPUT)  key  and  enter  the  date. 

6.  With  load  cell  assembled  and  unladen,  select  "kl"  (ADJ  B)  key. 
This  sets  the  bias  to  0.   To  check,  select  "k3"  (WEIGH).    If  a  weight 
within  .003  g  of  O.Og  is  not  obtained,  re-select  ADJ  B  until  satisfactory. 
Accurate  readings  should  also  be  verified  using  calibrated  weights. 

Filament  Loading 

1 .    Remove  Plexiglass  cover  from  sample  rack. 
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2.  Without  disturbing  adjacent  filaments,  lower  stabilizer  bolts  as  far  as 
possible  on  stations  to  be  loaded. 

3.  Adjust  the  spring  so  that  the  flag  will  operate  the  optical  trigger 
without  rubbing  against  it. 

4.  Remove  the  weight  from  the  station  to  be  loaded. 

5.  Using  tweezers,  grasp  the  desired  filament  by  one  end  of  it's  carrier 
and  hang  it  on  the  hook  attached  to  the  flag. 

6.  Using  tweezers,  hang  the  weight  on  the  bottom  hole  of  the  carrier. 

7.  Verify  the  flag  does  not  rub  against  the  trigger. 

8.  Lower  the  load  cell  elevator  far  enough  that  it  will  fit  onto  the  rack 
without  disturbing  the  hung  filament. 

9.  Without  anything  touching  the  top  of  the  load  cell,  press  "ADJ  B"  on 
the  HP-85. 

10.  When  the  screen  has  stopped  scrolling,  raise  the  elevator  until  the 
compliance  spring  is  completely  retracted  and  the  paper  carrier  is 
carrying  exactly  zero  load. 

11.  Verify  that  the  PC  records  the  station  as  unloaded  as  the  flag 
activates  the  optical  trigger. 

12.  Select  "LOAD"  on  the  HP-85  and  answer  the  cues: 
"ENTER  THE  SAMPLE  #" 

"ENTER  THE  STATION  #" 

13.  Select  "CONT" 

14.  Verify  the  weight  of  the  vial  shown  by  the  computer  is  correct  for 
the  station  being  loaded. 

15.  Ensure  the  direction  switch  on  the  elevator  is  in  the  "down" 
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position  and  the  power  switch  "off. 

16.  Check  the  setting  on  the  hot  wire  power  supply  and  turn  it  on. 
The  wire  should  glow  a  dull  red,  not  bright. 

17.  Without  touching  the  filament  or  disturbing  adjacent  filaments, 
carefully  burn  through  both  sides  of  the  paper  carrier  at  the  bottom  of 
the  slot. 

18.  Turn  off  the  hot  wire  and  return  it  to  it's  holder. 

19.  Select  "CONT"  on  the  HP-85  and  move  the  elevator  power 
switch  to  "on".   Monitor  the  elevator  as  it  lowers  to  ensure  it  does  not 
bottom  out  prior  to  being  unloaded. 

20.  Verify  the  PC  records  the  station  being  loaded  as  the  flag  blocks 
the  optical  trigger. 

NOTE.   If  the  filament  breaks  before  it  is  completely  loaded,  stop 
the  elevator,  record  the  mass  displayed  by  the  HP-85,  and  start  again 
from  step  24. 

21.  When  the  weight  hangs  free,  the  HP-85  will  beep  twice.   When 
this  happens,  stop  the  elevator.  After  two  minutes,  the  HP-85  will  beep 
twice  again,  allowing  you  to  continue. 

22.  Lower  the  elevator  until  it  can  be  withdrawn  without  disturbing  the 
filament. 

23.  Once  the  elevator  has  been  withdrawn,  carefully,  without 
disturbing  the  filament,  adjust  the  stabilizer  screw  so  that  it  is 
centered  on  the  vial  without  touching  it  and  has  a  vial/stabilizer  gap  of 
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not  more  than  1  mm. 

24.  Load  a  sheet  of  paper  into  the  plotter  and  select  "shift  k2"  (PLOT) 
on  the  HP-85,  answering  the  cues  as  required. 

25.  Record  filament  load  information  in  Fiber  Life  Data  Book  and  file 
load  plot  in  the  plot  binder. 

26.  Repeat  from  step  2  as  necessary. 

27.  When  all  desired  stations  have  been  loaded,  verify  all  stabilizer 
bolts  are  positioned  correctly  and  replace  Plexiglass  cover. 

28.  Verify  PC  is  reporting  the  correct  status  for  the  stations  loaded. 
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APPENDIX   B 

EQUIPMENT    DESCRIPTION 

Data  for  the  experiment  was  collected  by  suspending  fiber  samples  on 
the  individual  stations  (512  total)  of  a  large  sample  board.   Load  and  fail 
status  was  actively  monitored  by  individual  infrared  switchesthrough 
ribbon  wire  which  fed  the  data  to  an  interface  box  communicating  with   a 
PC- AT.    (Figure  B-l). 

The  sample  board  was  made  up  of  512  stations  arranged  in  8  rows 
and  64  columns.   Each  station  consisted  of  a  mounting  point,  a  high 
compliance  spring  (k=0.635  mm/g),  an  optical  switch  and  trigger,  a 
filament  sample,   a  weight,  and  an  adjustable  stabilizer  screw.     The  board 
was  constructed  of  aluminum  channel  stock  and  fiberboard.   When  not 
being  loaded,  the  stations  were  covered  with  a  protective  sheet  of  1/4" 
PMMA. 

The  trigger  was  a  wire  paper  clip  with  the  inner  loop  removed.   The 
top  half  of  this  outer  loop  was  then  wrapped  in  black  vinyl  electricians 
tape.   The  fiber  sample,  mounted  in  its  paper  carrier,  was  hung  on  this 
trigger  which  in  turn  hung  from  the, bottom  of  the  spring.   With  the  weight 
attached  to  the  bottom  of  the  fiber  carrier,  the  black  tape  blocked  the 
optical  switch.   When  the  filament  broke,  the  weight  dropped  away  while 
the  spring  retracted  the  trigger  allowing  light  to  pass  through  the  paper 


28 


clip  thus  closing  the  optical  switch,  which  was  detected  and  recorded  by 
the  PC. 

The  weights  were  plastic  vials  filled  with  lead  shot  and  were  loaded  to 
3  weight  levels:  12.20,  13.15,  and  14.17  g.  Two  rows  of  AS-4  graphite, 
batch  no.  008,  and  two  rows  of  AS-4  batch  no.  019  were  loaded  at  14.17 
g.,  one  row  of  008  and  one  row  of  019  were  loaded  at  13.15,  and  one  row 
of  008  and  one  row  of  019  were  loaded  at  12.20  g.   It  was  decided  to 
double  the  number  of  high  weighted  samples  in  order  to  obtain  more  data 
points  prior  to  censoring  the  experiment. 

Filaments  were  permanently  mounted  in  paper  carriers  to  facilitate  a 
reliable  and  convenient  method  of  storage  and  handling.    Once  assembled 
in  it's  proper  station,  the  paper  carrier  was  severed  using  a  hot  wire 
instead  of  scissors  to  reduce  the  chance  of  sending  a  shock  through  the 
filament.  [Ref.  3] 

Each  of  the  512  optical  switches  was  connected  to  a  custom  designed 
interface  box.   This  box  contained  a  power  supply  for  the  optical  switches 
as  well  as  a  multiplexing  circuit  which  fed  status  information  to  a 
MetraByte  I/O  card  in  the  PC.   The  board  was  divided  into  64   8-bit 
words  and  software  was  written  to  cycle  through  each  of  the  words  and 
to  report  any  changes  in  the  value  of  the  word.   For  example,  if  stations  1- 
8  are  all  loaded  the  PC  reads  a  value  of  0000  0000.   If  station  2  fails,  the 
next  time  this  byte  is  read,  it  would  show  a  value  of  0000  0010.   The  PC 
would  then  record  the  failure  and  continue  monitoring  the  board.    A  back- 
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up  power  supply  was  connected  to  the  interface  box  and  PC  in  case  of  a 
filament  failure  during  a  power  outage.   PC  software  is  reproduced  in 
Appendix  D. 

Start  times  for  each  filament  were  also  recorded  by  the  PC  when  a  bit 
change  from  1  to  0  was  detected.    However,  consistently  accurate  start 
times  were  highly  dependent  on  a  function  of  individual  spring  compliance, 
precise  trigger  construction  ,  loader  motor  speed,  and  masses  of  the 
weights.   In  order  to  eliminate  this  variability,   the  loader  software  used  to 
program  the  HP-85  also  incorporated  a  timing  mechanism.   As  the  load 
cell  was  lowered,  the  HP-85  monitored  the  weight  on  the  cell.   When  zero 
weight  was  detected,  a  beep  sounded  and  a  two  minute  timer  was 
started.   The  load  cell  was  left  under  the  sample  until  expiration  of  the  two 
minute  time  period.   Should  a  sample  have  failed  during  the  first  two 
minutes  of  life,  its  lifetime  would  be  recorded  to  the  nearest  1/10  second. 
After  two  minutes,  such  precise  record  keeping  was  not  necessary  and 
timing  was  turned  over  to  the  PC. 

The  HP-85  was  programmed  to  read  the  load  cell  approximately  5 
times  per  second  during  the  loading  sequence  and  could  therefore  provide 
accurate  data  on  the  loads  as  they  were  applied.    After  loading,  this 
information  was  sent  to  a  plotter  and  a  record  of  each  load  made  for 
analysis.   If  a  filament  ruptured  any  time  during  the  load  sequence,  up  to 
the  expiration  of  the  two  minute  time  period,  the  HP-85  would  record  the 
load  on  the  filament  at  the  time  of  failure  or,  if  already  loaded,  the  life  of 
the  filament  up  to  the  first  two  minutes. 
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Filaments  were  loaded  by  using  a  150  g.  capacity  Sensotec  load  cell 
mounted  to  an  electrically  controlled,  hydraulically  actuated  elevator.   The 
electric  motor,  Western  Gear  model  P5B24R3  incorporated  a  built  in 
reduction  gear  and  was  run  at  7.5  volts  giving  a  start  to  finish  load  run  of 
about  30  seconds  yielding  170-220  data  points  as  the  load  on  the  filament 
increased.    The  motor  was  mechanically  connected  to  a  disposable  5  cc 
hypodermic  syringe  which  acted  as  the  master  cylinder.    The  slave 
cylinder  was  an  identical  syringe  with  the  load  cell  mounted  on  a  platform 
secured  to  the  plunger.    This  arrangement  was  in  turn  mounted  to  a 
Plexiglass  platform  which  provided  a  stable  base  during  loading.   The 
slave  and  master  were  connected  with  3  feet  of  1/8"  ID  vinyl  tubing;  the 
hydraulic  fluid  was  tap  water. 
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AS4-019:  14.17  g. 
AS4-008:  14.17  g. 


AS4-019:  12.20  g. 
AS4-008:  12.20  g. 


AS4-019:  13.15  g. 
AS4-008:  13.15  g. 


AS4-019:  14.17  g. 
AS4-008:  14.17  g. 


Interface 
Box 


PC 


Figure  B-l.   Diagram  of  Fiber  Life  Testing  Apparatus. 
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APPENDIX    C 
ULTRAVIOLET   (UV)   LIGHT  SHIELDING 

The  life  of  certain  fibers,  (Kevlar  for  example)  are  known  to  be 
affected  by  exposure  to  UV  radiation.  To  assure  these  life  test 
experiments  were  performed  in  an  inert  environment,  an  investigation  into 
the  presence  of  UV  light  was  conducted. 

The  sample  rack  was  assembled  in  a  secure  laboratory  environment 
on  a  basement  level  with  windows  facing  Northeast.   The  rack  itself  was 
mounted  on  the  Southeast  wall  and  therefore,  could  not  receive  any  direct 
sunlight.    Windows  were  covered  by  opaque  commercial  grade  drapery. 
Room  lighting  was  fluorescent  overhead  with  plastic  covers.   The  sample 
racks  were  covered  with  1/4  inch  Poly  methyl  methacrylate  (PMMA) 
Plexiglas.   Light  from  the  two  sources  was  filtered  twice  before  reaching 
the  samples.   The  indirect  sunlight  passed  through  standard  grade  soda 
lime  window  glass  and  then  through  the  PMMA  rack  covers.    Overhead 
lighting  was  filtered  through  a  plastic  optical  diffuser  and  then  through  the 
PMMA  rack  cover. 

Tests  for  the  intensity  of  300-400  nm  UV  light  were  conducted  using  a 
100  Watt  mercury  lamp  at  a  distance  of  30  ft.  from  the  spectrometer.   All 
measurements  were  made  at  the  surface  of  the  sample  material.    Results 
of  the  spectrum  analysis  are  shown  in  Figures  C-l  through  C-3.   To 
explore  materials  for  improved  UV  shielding,  analysis  was  also  done  on 
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Polycarbonate  (Lexan)  and  a  UV  filtering  grade  of  PMMA  (Acrylite  OP- 
2).   Results  from  these  two  materials  are  shown  in  Figures  C-4  and  C-5. 
The  improvement  in  shielding  appeared  to  be  a  small  range  of  an 
additional  25  nm  to  that  already  provided  by  the  PMMA  and  was 
therefore  not  deemed  to  be  significant. 

Based  on  these  measurements,  the  amount  of  UV  light  reaching  the 
actual  fiber  samples  was  judged  to  be  insufficient  to  affect  life  testing. 
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Figure  C-l.   Results  of  UV  Transitu ttance  test  for  Soda  Lime  Glass. 
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Figure  C-2.   Results  of  UV  Transmittance  test  for  PMMA. 
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Figure  C-3.    Results  of  UV  Transmittance  test  for  Installed  Flourescent 
Lighting  Cover. 
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Figure  C-4.    Results  of  UV  Transmittance  test  for  Polycarbonate. 
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Figure  C-5. 


Results  of  UV  Transmittance  test  for  UV  Filtering  PMMA. 
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APPENDIX    D 


MAIN  PROGRAM: 


•**  Filamod7  ** 

1  Jim  Nageotte  Nov  13,  1990 

Start: 

DIM  newrd% (64) ,  OldRd%(64),  row%(8), 

Status$(2),  StarTime$(512) 

row%(l)  =  1:  row%(2)  =  2:  row%(3)  =  4: 

row%(4)  =  8:  row%(5)  =  16:  row%(6)  =  32 

row%(7)  =  64:  row%(8)  =  128 

a%  =  0 

OUT  &H303,  &H10  ' 

TYPE  Status2 

OldRd  AS  INTEGER 
END  TYPE 
DIM  Stat  AS  Status2 


*  Dimension  Variables 

*  Declare  constants 


*  Reset  I/O  Card 

*  Declare  file  type 

*  for  status  file 


TYPE  Times 

StarTime  AS  STRING  *  16 
END  TYPE 
DIM  Tim  AS  Times 


*  Declare  file  type 

*  for  Start  times  file 


CLS 

PRINT  "Ver  11.13.90" 

GOSUB  Gettime 

PRINT  "Power  up  at  "  ; 

DATE$ 

PrintStat  =  INP(&H3BD) 


TIME$;  "  on 


IF  PrintStat  =  223  THEN 
LPRINT  "Power  up  at  " 
"  on  " ;  DATE$ 
LPRINT  :  LPRINT 

END  IF 


TIME$; 


*  Print  start  up  message 

*  &  current  date  on  screen. 

*  Read  printer  status. 

*  223  indicates  that  the 

*  printer  is  ready. 

*  Print  start  up  message 

*  only  if  printer  is  ready. 


OPEN  "Status. Dat"  FOR  RANDOM 

AS  #1  LEN  =  2       ' 

FOR  idx  =  1  TO  64  ' 
GET  #1,  idx,  Stat 
01dRd%(idx)  -  Stat. OldRd 

NEXT  idx 

CLOSE  #1 

OPEN  "Times. Dat"  FOR  RANDOM  AS  #2 
LEN  =  16       ' 
FOR  i  =  1  TO  512 

GET  #2,  i,  Tim 

StarTime$(i)  =  Tim. StarTime 
NEXT  i 


*  Open  Status  File  for 

*  random  access. 

*  Read  previous  status 

*  from  disk  file. 


*  Close  Status  File. 

*  Open  Start  Times  File 

*  and  read  start  times 

*  for  all  512  stations. 


40 


CLOSE  #2 


*  Close  Start  times  file. 


DO 


GOSUB  Checksw 
FOR  i  =  1  TO  64 

IF  newrd%(i)  =  OldRd%(i)  THEN  GOTO  Same 
a0%  =  01dRd%(i)  :  al%  =  newrd%(i) 
FOR  j  =  1  TO  8:  jdx%  =  9  -  j 
b0%  =  FIX(aO%  /  row%(jdx%)) 
bl%  =  FIX(al%  /  row%(jdx%)) 
IF  b0%  =  bl%  THEN  GOTO  Nextbit 
Station  =  (i  -  1)  *  8  +  jdx% 
IF  bl%  <  0  OR  bl%  >  1  THEN 

PRINT  "Error  in  Status  bl%" 
BEEP 
END  IF 

IF  bl%  =  1  THEN 
PRINT  "Station  #";  Station; 
"Loaded   " ;  RdTime$;  "  on  "; 
RdDate$ 
GOSUB  Store 

PrintStat  =  INP(&H3BD) 
IF  PrintStat  -  223  THEN 

LPRINT  "Station  #"; 

Station;  "    Loaded  "; 

LPRINT  RdTime$;  "  on  "; 

RdDate$ 
END  IF 
BEEP 
END  IF 


*  Read  8  bit  status  words, 

*  until  something  changes, 


*  Compares  Status  read  from 

*  disk  to  that  read  from 

*  the  switches.  If  they  are 

*  not  equal  it  looks  more 

*  closely  at  the  switches. 

*  0  or  1  are  the  only  valid 

*  values  for  bl%. 


*"1"  indicates  "Loaded" 

*  Print  loaded  message  to 

*  screen. 

*  Store  new  Start  time  to 

*  disk  file. 

*  Check  printer  Status. 

*  If  printer  is  ready 

*  then  print  loaded 

*  message  to  printer. 


IF  bl%  -  0  THEN 
GOSUB  Elapsed 

PRINT:  PRINT  "Station  #••;  Station; 
"   ";  StarTime$ (Station) ; "to" ;PRINT 
LEFT$(RdTime$,  5);  "  " ;RdDate$: PRINT 
"  ";  YEARS;  "Years,  ";  DayT;  "Days  " 
PRINT  Min;  "  Minutes":  PRINT 
PrintStat  =  INP(&H3BD) 
IF  PrintStat  =  223  THEN 

LPRINT  "  ":  LPRINT  "Station  #"; 

Station;"   "; 

StarTime$ (Station) ;"  to  "; 

LEFT$(RdTime$,  5);"  ";  RdDate$ 

LPRINT  "  "; YEARS;  "Years,  ";  DayT;  "Days 

LPRINT  Min;  "Minutes" 

LPRINT  "  " 
END  IF 


*  "0"  indicates  failed. 

*  Find  elapsed  time. 

*  Print  failure  information 

*  to  screen. 


*  Check  printer  Status. 

*  If  printer  is  ready, 

*  then  print  failure 

*  information  to  printer. 
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END  IF 


Nextbit: 

aO%  =  aO%  -  bO%  *  row%(jdx%) 
NEXT'  j 


al%  =  al%  -  bl%  *  row%(jdx%) 


Same: 


01dRd%(i)  =  newrd%(i) 

OPEN  "Status.dat"  FOR  RANDOM  AS 

#1  LEN  =  2 

Stat.OldRd  =  01dRd%(i) 

PUT  #1,  i,  Stat 

CLOSE 

NEXT  i 


*  Update  disk  Status  file, 


LOOP 
END 


*  End  of  loop  that  compares 

*  old  status  to  new  and  end 

*  of  main  program. 


SUBROUTINES: 


Check  Switch  compares  the  current  status  to  the  status  read  off  the  disk 
file.  It  loops  until  there  is  a  change.  It  then  sets  a  flag  and 
returns  with  the  current  time  and  date  to  start  checking  each  of  the  512 
stations  and  determine  what  type  of  change  has  taken  place.  By  checking 
the  status  in  8  bit  words  and  in  such  a  small  loop,  we  keep  the  time 
required  to  check  all  switches  and  the  response  time  to  a  minimum. 


Checksw: 


DO 

GOSUB  Gettime 

a%  =  0:  idx  =  0:  RdTime$  =  TIME$ 
RdDate$  =  DATE$ 

FOR  i  =  1  TO  8:  OUT  &H301,  row%(i) 
FOR  j  =  0  TO  7:  OUT  &H302,  j 
idx  =  idx  +  1 
newrdl(idx)  =  INP(&H300) 

IF  newrd%(idx)  <>  OldRd%(idx) 
THEN  a%  =  1 
NEXT  j 
NEXT  i 

LOOP  UNTIL  a%  =  1   ' 
BEEP 
RETURN 


*  Get  current  time. 

*  Get  current  date. 

*  Select  Row  to  read. 

*  Select  word  to  read. 


*  Set  Flag  if  Status  has 

*  changed. 


*  Loop  until  Status  has 

*  changed. 
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Elapsed  Time  calculates  the  time  between  the  Start  time  read  from  disk 
and  the  current  time.  It  calls  two  other  sub  routines  to-do  this,  Julian 
and  Minutes.  Julian  is  called  twice,  the  first  time  it  converts  the 
Start  time  to  a  Julian  date  and  the  second  time  it  converts  the  current 
date  to  a  Julian  date.  The  conversion  is  necessary  to  allow  subtraction 
to  be  performed  on  the  dates.  Having  arrived  at  the  number  of  Years  and 
Days  before  the  fiber  ruptured,  we  call  the  subroutine  Minutes  to 
subtract  the  times  involved  and  return  with  the  correct  number  of 
minutes.  An  adjustment  to  the  number  of  days  may  be  made  based  on  the 
number  of  minutes. 


Elapsed: 

ConvDate$=RIGHT$ (StarTime$ (Station) 

,10) 

GOSUB  Julian 

StartYr  =  Year 

StartDy  =  Days 

ConvDate$  =  RdDate$ 

GOSUB  Julian 

StopYr  =  Year 

StopDy  =  Days 

YEARS  =  StopYr  - 

StartYr 

DayT  =  StopDy  - 

StartDy 

IF  DayT  <  0  THEN 

DayT  =  DayT 

+  365 

YEARS  =  YEARS  -  1 

END  IF 

GOSUB  Minutes 

IF  Min  <  0  THEN 

Min  =  Min  + 

1440 

DayT  =  DayT 

-  1 

END  IF 

RETURN 

Julian: 

Mon  =  VAL(LEFT$(ConvDat 

e$,  2)) 

Day  =  VAL(MID$(ConvDate$,  4,  2)) 

Year  =  VAL(RIGHT$ (ConvDate$,  2)) 

SELECT  CASE  Mon 

CASE  1 

Days  = 

Day 

CASE  2 

Days  = 

Day  + 

31 

CASE  3 

Days  = 

Day  + 

59 

CASE  4 

Days  = 

Day  + 

90 

CASE  5 

Days  = 

Day  4- 

120 

CASE  6 

Days  = 

Day  + 

151 

CASE  7 

Days  = 

Day  + 

181 

*  Select  Start  Date. 

*  Convert  to  Julian  date, 

*  Select  Current  Date. 

*  Convert  to  Julian  date. 


*  Subtract  years. 

*  Subtract  days. 


*  Calculate  minutes. 


*  Get  month,  day  and  year 

*  from  the  date  string. 


*  Calculate  total  days. 
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CASE  8 

Days  =  Day  +212 
CASE  9 

Days  =  Day  +  24  3 
CASE  10 

Days  =  Day  +  27  3 
CASE  11 

Days  =  Day  +304 
CASE  ELSE 

Days  =  Day  +  3  34 
END  SELECT 
RETURN 

Store  simply  opens  the  Start  Time  file  on  disk,  a  random  access  file, 
and  stores  only  the  new  time  without  reading  or  re-writing  the 
unaffected  times.  The  file  is  then  closed  and  the  routine  exited. 

Store; 

OPEN  "Times. Dat"  FOR  RANDOM  AS  #2  LEN  =  16 

Tim.StarTime  =  LEFT$ (RdTime$,  5)  +  "  "  +  RdDate$ 

PUT  #2,  Station,  Tim 

CLOSE  #2 

StarTime$ (Station)  =  Tim.StarTime 

RETURN 

Minutes  recovers  the  values  of  Minutes  and  Hours  from  the  time  strings 
and  converts  them  to  Minutes  only  for  ease  of  mathematical  operations 
in  the  time  keeping. 

Minutes: 

HrSt$  =  LEFT$ (StarTime$ (Station) ,  2) 

Min$  =  MID$ (StarTime$ (Station) ,  4,  2) 

Hrs  =  VAL(HrSt$) 

MinSt  =  VAL(Min$) 

MinS  =  Hrs  *  60  +  MinSt 

Hrf$  =  LEFT$(RdTime$,  2) 

MinF$  =  MID$(RdTime$,  4,  2) 

Hrf  =  VAL(Hrf$) 

MinFn  =  VAL(MinF$) 

MinF  =  Hrf  *  60  +  MinFn 

Min  =  MinF  -  MinS 

RETURN 

Gettime  reads  I/O  ports  to  read  the  time  directly  from  the  system  clock. 
It  then  converts  these  values  to  a  string  that  resembles  the  "TIME$" 
function.  The  program  then  uses  this  string  called  '^tin^"  for  Real 
Time  String. 
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Gettime: 

Sec$  =  LTRIM(STR$(INP(&H2C0)  AND  &HF) ) 

Secs$  =  LTRIM(STR$(INP(&H2C1)  AND  &HF) ) 

Minute$  =  LTRIM(STR$ (INP(&H2C2)  AND  &HF) ) 

Minutes$  =  LTRIM(STR$ (INP(&H2C3)  AND  &HF) ) 

Hour$  =  LTRIM(STR$(INP(&H2C4)  AND  &HF) ) 

Hours$  =  LTRIM(STR$(INP(&H2C5)  AND  &HF) ) 

Rtime$  =  "00: 00; 00" 

MID$(Rtime$,  1)  =  Hours$ 

MID$(Rtime$,  2)  =  Hour$ 

MID$(Rtime$,  4)  =  Minutes$ 

MID$(Rtime$,  5)  =  Minute$ 

MID$(Rtime$,  7)  =  Secs$ 

MID$(Rtime$,  8)  =  Sec$ 

RETURN 
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Read  Data  Port 


Determine  which 
station  has  changed 


If  loaded,  display 
Sta.  No.  and  start  time/ date. 
Record  start  time/  date 
on  floppy  drive. 


If  not  loaded,  display  Sta.  No. 
and  start/  fail  time  and  date. 
Record  Sta.  No.  and  fail 
time/  date  on  floppy  drive. 


Figure  D-l.   Flowchart  for  Life  Test  Monitoring  Program. 
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APPENDIX    E 


I.     LIKELIHOOD  DERIVATION  FOR  EXACT  DATA 


For  Exact  Data: 


Lexact  =  f(xr  a>  P)  =  <^dx  F(xi'  a'  P) 


:F(x,a,p)=l-exp|-^a 


Taking  the  derivative  with  respect  to  x: 


f(x)=f=-e«i 


OC] 


-rrax^ 


at  x 


a-1 


PvP 


enix  -  -9- 


a] 
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The  Likelihood  for  all  data,  x;,  = 


L(xi)=n)f(xi,a,P) 


n 

n 
i=i 


.a/ x. 

pip 


a-1 


e«p 


(  *•  A 


ai 


(1) 


Rearranging, 


n 

n 
i=i 


-1  (1-cc) 
ap    p  e«p 


(X.  > 

i  P\ 


ai 


n 

n 
i=i 


a 


-a     (a-1) 
ap      x.  e«p| 


(  x.  >\ 
1 


V  P  )    J 


nD-na  "       (a-1) 
a   p  U  x.  e«p 


ex.  \ 
i 

P  J 


ai 


nQ-na  "       (a-1)         f       -a  "     a] 
a   p  n  x.  enpJ  -  p       I  xu 

i=1    !  1  i=1    l\ 


(2) 
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Let  A  =  InL  =  In  f(Xj,a,p) 


=  ln 


a"p-n"flxi(a-,)eHp{-p-aix1;) 


n  Ina  -  na  lnp  +  Xln 


:°-»eKP{-rix°}] 


n  11        r  f~  -     "  ~1 

=  n  Ina  -nalnp  +  (a-  DlnXx  +  XlrJenp-J  -  p    Ix°  > 

1=1         1=1    L  L  ,=i     J 


n  _a  n 

n  lna-nalnB  + (a-  1)ln  I  x.  -  B       Ix 

.=1    »  i=1    ' 


(3) 


To  find  the  maximum  Likelihood,  (Lmax  ),  set  the  derivative  of  A  with 
respect  to  a  and  p,  to  zero: 


3A  _  3A  =0 

3a  "'  ap 
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Taking  the  partial  derivative  of  eq.  (3)  with  respect  to  a  and  p  and 
setting  it  to  zero  yields: 


^-  =  -£--nlnp+  Ilnx.  +  (3  aln(3Xxa-p  aIxalnx   =0 
ca      a  !  j  i=1  1        i 


(4) 


3A  __  na      __a_y„a  _  n 

ap  "     P     p(u+,)..  ' 


(5) 


Rearranging  eq.  (5): 


n  a  a 


>,?ix' 


1    n 


(6) 


i    n       "\ 


1=1    v 


(7) 
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Now  replace  the  n  in  the  2nd  term  of  eq.(4),  with  eq.  (6): 


.-a 


.-a 


77-f4r£xaW+  Ilnx.  h[T    nplxa-p  "Xxalnx  =0 
a      U   i=1    'J  i=1  i=1   «  1=1   ' 


and  cancel  with  the  4th  term: 


77+  Ilnx  -p"aIxalnx  =  0 
a     .=  i       '  i=i  ' 

From  eq.(6),  substitute 


n(3    =  Ixa 

and  divide  both  sides  by  n  yielding: 


1       1  n 
a+J-IInx  - 

i  =1 

Ixalnx 

.=i  ' 

n 

Ixa 

,=  1     ' 

0 
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Rearranging  to  set  the  right  side  equal  to  a: 


Ixalnx 

I  .    .     l  i  *n 

,    1 

i  =1  ,or 


a  = 


n 

Zxalnx 

,=  i   ' 

1      " 

--J-  Ilnx 

n,=i       ' 

n 

Ixa 

i=1       ' 
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(8) 


APPENDIX    F 


I.     LIKELIHOOD  DERIVATION  FOR  EXACT  DATA  WITH 
TYPE  I  AND  TYPE  IV  CENSORING 


For  Exact  data  which  has  both  Type  I  and  Type  I-v  (Variable:  Cause  of 
failure  not  quantifiable  in  this  experiment)  censoring: 

Likelihood,  =  [nffx.'^P)]!!!1  "  F(\-a,p)][l  -  F(xr  a-P)]  (n"m) 
Exact  Type  I-v  Type  I 


ffcAfl^l-eJ-fe 


l-F(x\,a,P)  =ex  J- 


l-F(Xj,a,P)=exri-^ 


A 

:        n    =The  total  number  of  fibers 

m   =The  total  number  of  fibers  that  failed 
m    =The  numberof  fibers  failing  in  life 

m  -  m    =The  numberof  fibers  subject  to  Typel  -  v 
censoring 
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Therefore, 


L 


_  m,  n  -m,a  I    I      a-l„  ft-a     ct  1 

a  '(3       }|x     exp{-p    x,  } 


riexp  - 


'O" 


vp; 


expi- 


(n-m) 


(D 


Let  A  =  InL 


m,ami 


=  ln  a    '(3     '   nx^expl-  (TV} 

i=1 


in 


f    X 


ln|   n  e«p 

i  =m  . 


+  In 


(n-m)" 


expi- 


m  m 

1  -a     ' 


=  m,lna -m,a  lnS  + (a-  1)  Ilnx   -  B     Ixa 
i  i        r  i       '  i    ' 

i=i  i=i 


P       X  x     -  (n  -  m) 


(  x.  \ 
j 

(3  J 


(2) 
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To  find  the  maximum  likelihood  of  L, 

Set^=^  =  0 
3a      dp 

Taking  the  partial  derivative  of  eq.(2)  with  respect  to  a  yields: 


dA       mi 


m 


^  =  -F7i-m,lnp+Ilnx-p     Xxalnx  -{3     lnplxa 
3a       a  '  i=i       '  1      i=i   '        '  i=i   V 


m 


rv      m  n  r,  m  n 

I  x.  Inx.-B     InB  I  x 

i        i  .         i 


i  =m 


i=  m 


X.A 


-(n  -m)M 


In 


^       =0 


(3) 


Taking  the  partial  derivative  of  eq.(2)  with  respect  to  p  yields: 


dA 


m 


in 


m,aP"1  +  a(3"a"1  Xxa  +  a(Ta_1  I  x  "  +  aB""_1(n  -  m)xa  =  0 

i  =1  i  =m  , 


a+1 


Multiplying  through  by  , 


m 


m 


-mBa+Ixa+    I  £a  +  (n  -m)xa  =  0 
1  ,=i    '       i-m.    '  J 
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Now  divide  by  rrij  and  rearrange  to  isolate  (3: 


K        m 


rm 


l  m 


Ixa+    X  x    +  (n  -  m)xa 

j 


i  =1  i  =m  . 


p= 

1 

ml 

rm'                m        a                                T 

Ixa+    X  x     +(n  -m)xa 

.   ,    i       .          i                        j 

;i=1          1=111,                             J_ 

i 

"a 

(4) 

To  solve  for  a,  first  multiply  eq.  (3)  by  (3a  and  rearrange: 


o/m,  m'  1  "'  m        a 

P    -sL-m1lnp+  Ilnx.   +lnpix«  +  lnp  I  x!  +(n  -m)xalnp 

V  1=1  /  1=1  i=m  . 


in 


m 


x-^      ex.  ■  x-"     a  u  i      a  \      ix  t 

L  x   Inx    +    X  x    Inx    +  ( n  -  m ) x    I n x 

,ii  ii  ij 

i=l  i=  m  .  J 


Now  gather  the  lnp  terms. 


m 


m. 


P     -F7i-m,lnp+  Ilnx 


i  m 


+  lnfl  Ixa+    I  x.  +(n  -m)xa 


i=1  i=m 


in 


L  x   Inx    +    L  x    Inx    +(n-m)x   Inx 


,ii  ii 

i=l  i=m . 


in 
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Divide  both  sides  by  mj  and  rearrange  to  bring  a  to  the  left  hand  side: 


1  1 

a  n,      ft 


m 


m 


£xalnx.+    X  x  .  Inx  +  (n  -  m)xalnx 
.  ,   i        i  i         i  j        j 

i=1  i=m.  J 


m 


m 


Ixa+    I  £u  +  (n-m)xa 

i=1  i=m.  J 


i.     + 


Inp 


i      ' 
lnB  +  7^-1  Inx. 

mi,=  i       ' 


(5) 


Obtain  from  eq.  (4), 


m 


m 


m,6    =  Ixu+    X  x    +  (n  -m)xa 

i  =1  i  =m  , 


(6) 


and  substitute  for  (a)  in  eq.  (5): 


1  1 

a  n ' 


in 


in 


i=1  i  =m 


a.      ^ 

Inx 

i         i 


Ixalnx .+    I  x    Inx   +(n-m)xalnx 


in 


[mp*)np 


J    S>  + 


i      > 
lnB  +  — —  Ilnx 

mi,=  i       ' 
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Cancel  the  ln(3  terms  and  invert  both  sides: 


a  = 


m 


i=1 


1  =m 


Sxalnx .+    Z  x.  Inx.  +  (n -m)xalnx 


j        j 


m,P 


1    m' 
+  — Xlnx 

mn=i       » 


Finally,  eliminate  (3  by  substituting  eq.  (6)  for  the  denominator: 


a  = 


i  m 

xlnx   +    2,  x    Inx   +  ( n  -  m )  x  .  I  n  x 

1  j        j 


i=l  i  =m 


m 


m 


x.   +  (n  -  m)x     +    L  x  . 

1    '  J  ' 

i  =1  J        1=  m 


1     m 

mn=i       ' 
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Figure  1 .  Relationship  of  Fiber  sample  sets  from  Spool  AS4-008 
population. 
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Figure  2.  Relationship  of  Fiber  sample  sets  from  Spool  AS4-019 
population. 
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Figure  3.   Expected  Rank  of  008  strength  set  {Ds} 
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Figure  4.   Expected  Rank  of  019  strength  set  {D5} 
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Figure  5.   Expected  Rank  of  008  strength  sets  {D5}  and  {Dts} 
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Figure  6.   Expected  Rank  of  019  strength  sets  {Ds}  and  {Dts} 
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A 

B 

C 

1 

AS4-008 

2 

3 

4 

5 

6 

xi  (Ds|Ns) 

xi    (Dts|Nts) 

F*({Ds+Dts}|N 

7 

5.909 

-6.277309565 

8 

6.393 

-5.382913936 

9 

6.666 

-4.917448539 

10 

6.711 

-4.600324093 

11 

7.12 

-4.359337455 

12 

7.997 

-4.164810785 

13 

8.017 

-4.001618953 

14 

8.54 

-3.860998686 

15 

8.9 

-3.737415596 

16 

9.085 

-3.627150043 

17 

9.364 

-3.52758132 

18 

9.58 

-3.436793307 

19 

10.037 

-3.35334254 

20 

10.04 

-3.276114499 

21 

10.047 

-3.204230672 

22 

10.136 

-3.136986293 

23 

10.234 

-3.073807343 

24 

10.263 

-3.014220084 

25 

10.551 

-2.957828979 

26 

10.564 

-2.904300374 

27 

10.588 

-2.853350225 

28 

1 0.592 

-2.804734725 
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-2.758243047 
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10.776 

-2.713691648 
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10.909 
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10.929 

-2.629785756 

33 

10.973 

-2.590164254 

34 

10.987 

-2.551943695 

35 

10.998 

-2.515024401 

36 

11.169 

-2.479316935 

37 

11.199 

-2.444740751 

38 

1 1 .279 

-2.41122305 

Figure  7a.   Expected  Rank  data  of  merged  008  strength  set  {D5+Dts} 
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A 
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39 

11.378 

-2.378697815 

40 

11.4 

-2.347104998 

41 

1 1 .41 7 

-2.316389819 

42 
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-2.286502166 
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Figure  7b.   Expected  Rank  data  of  merged  008  strength  set  {Ds+Dt5} 

(cont.) 
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Figure  7c.   Expected  Rank  data  of  merged  008  strength  set  {D5+Dt5} 

(cont.) 
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Figure  7d.   Expected  Rank  data  of  merged  008  strength  set  {Ds+Dts} 
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Figure  7e.   Expected  Rank  data  of  merged  008  strength  set  {D5+Dt5} 

(cont.) 
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Figure  7f.   Expected  Rank  data  of  merged  008  strength  set  {Ds+Dts} 

(cont.) 
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Figure  7g.   Expected  Rank  data  of  merged  008  strength  set  {D5+Dts} 

(cont.) 
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Figure  7h.   Expected  Rank  data  of  merged  008  strength  set  {Ds+Dts} 

(cont.) 
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Figure  7i.   Expected  Rank  data  of  merged  008  strength  set  {Ds+Dts} 

(cont.) 
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Figure  7j.   Expected  Rank  data  of  merged  008  strength  set  {Ds+Dts} 

(com.) 
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Figure  8a.   Expected  Rank  data  of  merged  019  strength  set  {Ds+Dts} 
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Figure  8b.   Expected  Rank  data  of  merged  019  strength  set  {Ds+Dts} 

(cont.) 
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Figure  8c.   Expected  Rank  data  of  merged  019  strength  set  {Ds+Dt5} 

(cont.) 


77 


E 

F 

G 

115 

14.01 

-1.127438617 

116 

14.069 

-1.116643872 

117 

14.085 

-1.105926941 

118 

14.16 

-1.095286311 

119 

14.213 

-1 .084720508 

120 

14.27 

-1 .074228096 

121 

14.271 

-1.063807675 

122 

14.305 

-1.053457881 

123 

14.328 

-1.043177384 

124 

14.332 

-1.032964885 

125 

14.346 

-1.022819119 

126 

14.353 

-1.012738851 

127 

14.37 

-1 .002722875 

128 

14.416 

-0.992770013 

129 

14.434 

-0.982879117 

130 

14.481 

-0.973049063 

131 

14.485 

-0.963278755 

132 

14.528 

-0.953567122 

133 

14.54 

-0.943913114 

134 

14.579 

-0.934315709 

135 

14.586 

-0.924773904 

136 

14.596 

-0.915286721 

137 

14.648 

-0.9058532 

138 

14.649 

-0.896472405 

139 

14.676 

-0.887143417 

140 

14.678 

-0.877865339 

141 

14.681 

-0.868637291 

142 

14.69 

-0.859458411 

143 

14.697 

-0.850327856 

144 

14.702 

-0.841244799 

145 

14.767 

-0.832208431 

146 

14.825 

-0.823217957 

147 

14.83 

-0.814272599 

148 

14.855 

-0.805371595 

149 

14.929 

-0.796514195 

150 

14.95 

-0.787699666 

151 

14.97 

-0.778927288 

152 

15.003 

-0.770196353 

Figure  8d.   Expected  Rank  data  of  merged  019  strength  set  {D5+Dts} 

(cont.) 
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Figure  8e.   Expected  Rank  data  of  merged  019  strength  set  {Ds+Dts} 

(cont.) 
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Figure  8f.   Expected  Rank  data  of  merged  019  strength  set  {D5+Dt5} 

(cont.) 
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Figure  8g.   Expected  Rank  data  of  merged  019  strength  set  {D5+Dt5} 

(com.) 
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Figure  8h.   Expected  Rank  data  of  merged  019  strength  set  {D5+Dts} 

(cont.) 
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Figure  8i.   Expected  Rank  data  of  merged  019  strength  set  {D5+Dts} 

(cont.) 
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Figure  8j.   Expected  Rank  data  of  merged  019  strength  set  {D5+Dts} 

(cont.) 
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Figure  8k.   Expected  Rank  data  of  merged  019  strength  set  {Ds+Dt5} 

(cont.) 
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Figure  10.   Likelihood  surface  and  contours  of  019  strength  set  {D5} 
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Figure  11.   5  %  and  1  %  confidence  intervals  for  008  strength  set  {Ds} 
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Figure  13.   5  %  and  1  %  confidence  intervals  for  008  strength  sets  {D5} 

and  {Dts} 
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and  {Dts} 
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Figure  15.   Expected  Rank  of  008  strength  sets  {D5  +Dts}  and  {D5} 
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Figure  16.   Expected  Rank  of  019  strength  sets  (D5  +Dts}  and  {D5} 


93 


punpLOOLopuopunpLnpunpuop 

r   ddoiaodedsN   c<dioid**eon 

C\ICMCMT-T-T-1-T-T-T-1-T-T-f-T-1-T- 

Figure  17.   5  %  and  1  %  confidence  intervals  for  008  strength  sets  {Ds} 

and  (Ds  +  Dts} 


94 


a  a  5 

CO.    CO.    "o 

J    J    CO 


I     i 


LO 
00 

o 

00 
LO 

r^ 
o 

LO 

CD 

o 

CD 

ID 
LO 

o 

LO 
LO 

o 

LO 
CO 

o 

CO 
LO 

c\i 

o 
c\i 

LO 


LO 

o 


pLOOLOOLOpLOp 

l-    ddoidodcoN    n 

OJCMOJt-t-t-t-t-i- 


LOpLOOLOOLOp 
CO      CO      LO      LO      ^ '^      CO      CO 


03 
.C 
Q. 

< 


ejag 
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Figure  20.   Expected  Rank  of  008  strength  sets  {Ds  +Dts}  and  {Dt2s} 
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Figure  21.   Expected  Rank  of  008  strength  sets  (Ds  +Dts}  and  {Dt3s} 
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Figure  22.   Expected  Rank  of  019  strength  sets  {D5  +Dt5}  and  {Dt]s} 
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Figure  23.   Expected  Rank  of  019  strength  sets  {Ds  +Dts)  and  {Dt2s} 
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Figure  24.   Expected  Rank  of  019  strength  sets  (Ds  +Dts}  and  {Dt3s} 
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Figure  25.   Expected  Rank  of  008  strength  sets  {Ds  +Dts}  and  {D} 
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Figure  26.   Expected  Rank  of  019  strength  sets  {Ds  +Dts}  and  {D} 
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Figure  27.   Likelihood  contours  of  008  strength  set  {D} 
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Figure  28.  Likelihood  contours  of  019  strength  set  {D} 
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Figure  29.  Likelihood  contours  of  008  life  set  {Dt]} 
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Figure  30.  Likelihood  contours  of  008  life  set 
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Figure  31.  Likelihood  contours  of  008  life  set  {Dt3} 
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